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ABSTRACT 

As part of a larger program aimed at better quantifying the uncertainties in stellar computations, we attempt to calibrate the extent 
of convective overshooting in low to intermediate mass stars by means of eclipsing binary systems. We model 12 such systems, with 
component masses between 1.3 and 6.2 M 0 , using the detailed binary stellar evolution code stars, producing grids of models in both 
metallicity and overshooting parameter. From these, we determine the best fit parameters for each of our systems. For three systems, 
none of our models produce a satisfactory fit. For the remaining systems, no single value for the convective overshooting parameter fits 
all the systems, but most of our systems can be well described with an overshooting parameter between 0.09 and 0.15, corresponding 
to an extension of the mixed region above the core of about 0.1-0.3 pressure scale heights. Of the nine systems where we are able to 
obtain a good fit, seven can be reasonably well fit with a single parameter of 0.15. We find no evidence for a trend of the extent of 
overshooting with either mass or metallicity, though the data set is of limited size. We repeat our calculations with a second evolution 
code, mesa, and we find general agreement between the two codes. For the extension of the mixed region above the convective core 
required by the mesa models is about 0.15-0.4 pressure scale heights. For the system El Cep, we find that mesa gives an overshooting 
region that is larger than the stars one by about 0.1 pressure scale heights for the primary, while for the secondary the difference is 
only 0.05 pressure scale heights. 
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1. Introduction 

Using observations to constrain stellar evolution models is 
one of the primary drivers of stellar astrophysics research. 
Nevertheless, it is not straightforward to directly link observa¬ 
tions and theory in a fully consistent manner, properly taking 
into account the uncertainties, particularly the theoretical ones. 
As a result, observers often consider their preferred set of stellar 
evolution tracks as intrinsically correct, and only recently com¬ 
parisons between the stellar parameters gathered from different 
sets of mo dels have begun to be per formed, mostly for low-mass 
stars (e.g. ICasagrande et al.ll20lH) . This “exercise” should be 
regularly performed, but in order to do that, particularly for a 
large number of stars, one needs to develop a dedicated auto¬ 
matic tool. 

To this end, we have set up the “Confronting Uncertainties 
in Stellar Physics’ (CUSP) project, the aim of which is to bet¬ 
ter quantify the impact of theoretical uncertainties in the use 
of stellar evolution models to determine fundamental stellar pa¬ 
rameters from observables. We intend to make use of bonnsaQ 
dSchneider et al.ll2014l) . a publicly available tool, which allows 
one to derive stellar parameters (e.g., mass, radius, age) from a 
set of observational parameters (e.g., effective temperature, sur¬ 
face gravity, rotational velocity), properly accounting for the ob¬ 
servational uncertainties. We will then explore the impact of both 
observational and theoretical uncertainties in the analysis of par¬ 
ticular astrophysical problems which are directly linked to the 
estimation of stellar fundamental parameters on the basis of evo¬ 


1 The bonnsai web-service is available at 

www.astro.uni-bonn.de/stars/bonnsai. 


lutionary tracks, such as the derivation of the masses and radii of 
transiting exoplanets. 


In its current status, bonnsai contains sets of evolutionary 
tracks only for intermediate to massive stars (M>5 M e ) calcu- 
lated for three different metallicities: Galactic, LMC and SMC 
(iBrott et al.ll201 It iKohler et al.ll2015T) . We aim at extending the 
model database to lower mass stars from 0.8 to 10 M 0 , and for a 
variety of different metallicities. This extension of bonnsai will 
be made by calculating n ew large grids of models initially with 
two d ifferent codes stars (Egg eton 


12009) and mesa dPaxton etall 


1971: IStancl iffe & Eldridgel 


2011), with the intention to add 


further grids calculated with other codes. However, the first step 
in the creation of these grids is the calibration of certain parame¬ 
ters. In this work, we focus on the issue of convective overshoot¬ 
ing. 


One of the key uncertainties in the evolution of main se¬ 
quence stars is the size of the convective core. It is widely ac¬ 
cepted that models based on either the Schwarzchild or Ledoux 
criterion produce cores that are too small to match observations 
and hence s ome form of oversh ooting must be applied to stellar 
models (e.g. lMaeder & Mevnetll 1 99~P) . By overshooting we sim¬ 
ply mean that the chemically mixed region in the star’s core has 
been extended beyond the convective boundary predicted from 
standard stellar theory. This additional mixing could be caused 
by any number of phenomena, and is not necessarily related to 
the motion of material driven by convection ‘overshooting’ the 
formal convective boundary. 


Convective overshooting is a key component of canonical 
stellar models but it requires calibration. Many possible methods 
for calibration exist, including the fitting of isochrones to stel¬ 
lar cluster colour-magnitude diagrams (e.g. VandenBerg et al.l 
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l200d ). More recently, asteroseismology has opened up a new 
aven ue for calibration of mixing properties in s t ellar interiors 
(e.g. lMontalban et~aPl2013l : lGuenther et al.ll2014h . lAertsl (120151) 
presents 16 OB dwarfs for which asteroseismic determinations 
of the extent of overshooting have been made. Similar anal- 
yses have also been c arried out by iNeiner et al.l (12012l) and 
iTkachenko et alJ (l2014lt . 

Here we focus only on calibrations using binary systems. 
[Schroder et al.l (119971 ) used / Aurigae-type systems to attempt 
to calibrate the extent of overshooting. They could find adequate 
fits to their systems (with masses between 2.5 and 6.5 M 0 ) us¬ 
ing overshooting equivalent to 0.24 - 0.32 pressure scale heights. 
Subseauentlv. lPols et al.l ( 1997 ) used the same evolutionary code 
and overshooting prescription to look at 49 eclipsing binary 
systems taken from the compilation of lAndersenl (1199 ll) . They 
found that models with and without overshooting could ad¬ 
equately fit the observations of the majority of the systems. 
However, for three systems (namely AI Hya, WX Cep and TZ 
For) models with enhanced mixing provided a better fit. 

Further attempts to calibrate convective overshooting have 
been made by Idaretl d2007l) . He used 13 double-lined eclips¬ 
ing binaries covering a range of evolutionary states and masses 
(from 1.3 to nearly 30 M 0 ). Moderate amounts of overshooting 
of around 0.2 times the pressure scale height were found to best 
fit the data, with little evidence for a mass de pen dency. This is 
in contrast to the earlier work of iRibas et all ( l2000h . who sug¬ 
gested a mass dependence for overshooting may exist, based on 
a samp le of 8 stars between 2 and 12 M Q . More recent work by 
iMeng & Zhand d2014l) . using four eclipsing binaries in the mass 
range 1.3-3.6M 0 , also finds no evidence for a mass dependency 
to overshooting, though the mass range is much smaller than the 
two studies mentioned above. 

In this work, we revisit the issue of calibrating convective 
overshooting using eclipsing binaries. Ultimately, the aim is to 
arrive at a reliable determination that can be used for the compu¬ 
tation of large grids of low-mass stellar models for use with the 
bonnsai tool. 


2. Stellar models 


Computations in this work were made using the stars 
stellar evolution code originally developed by Eggletor 
and updated by many authors (e.g. iPols et al 
The code is freely available for download from 
http://www.ast.cam.ac.uk/~stars. This code solves the 
equations of stellar structure and chemical evolution in a fully 
simultaneous manner, iterating on all variables at the same 
time in order to converge a model (see IStancliffel I2006L for 
a detailed discussion ). The version employed here is that of 
[Stancliffe & Eldridgel ( 2009 ) which was developed for doing 
binary stellar evolution. The code treats all forms of mixing by 
means of a diffusive formalism (lEggletonlll972l) . 

ISchroder. Pols. & Eggletonl d 19971) describe the implemen¬ 
tation of overshooting in the code. Rather than applying an ex¬ 
tension to the convective region that is some fraction of a pres¬ 
sure scale height, this implementation makes an adjustment di¬ 
rectly to the convective criterion (in this case the Schwarzchild 
criterion). A region is determined to be convectively unstable if 
Vrad > V ad - <5, where 


6 = 


_£ov_ 

2.5 +20/+ 16/ 2 


( 1 ) 


with / being the ratio of radiation to gas pressure and d ov is a 
constant that must be determined (i.e. the overshooting param¬ 
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Fig. 1 . Hertzsprung-Russell diagram showing eclipsing binaries 
from the [Torres et al. ( 201 0) sample, with the systems we have 
selected shown in red. Stellar evolutionary tracks computed 
without overshooting are shown to give an approximate indica¬ 
tion of mass and evolutionary state. 


eter). This formalism ensures a smooth transition between sys¬ 
tems with and without convective cores. 

Each model is evolved from the pre main sequence using 999 
mesh pointfl The mixing length a is set to 2.0, based on cali¬ 
bration to a Solar model. No mass loss is included. In addition, 
each binary is placed in a wide orbit so that there is no interaction 
between the components - we do not attempt to reproduce the 
observed orbital period. Similarly, these models are non-rotating 
and we do not try to reproduce the observed rotational velocity 
for those systems where it has been measured. Rotational mix¬ 
ing would also act to increase the size of the chemically mixed 
region. Provided our target systems do not rotate rapidly, the use 
of non-rotating models should suffice. Each system is evolved to 
the point where the primary has comfortably exceeded the ob¬ 
served primary radius. 


3. Results 

We attempted to model 11 binary systems: V364 Lac, AI Hya, 
El Cep, TZ For, WX Cep, VI031 Ori, SZ Cen, AY Cam, AQ 
Ser, V539 Ara and CV Vel. The properties of each of these sys¬ 
tems are listed in Table[3] These particular syst ems were selected 
from the sample of lTorres. Andersen. & Gimenezl ( 2010 ). which 
lists the known eclipsing binaries whose parameters have been 
determined to better than 3%. We chose those systems which ap¬ 
pear the most evolved in the Hertzsprung-Russell diagram (see 
Fig. |T}, with both components of the system being clearly sepa¬ 
rated from the zero age main sequence. Evolved systems should 
be more sensitive to the effects of overshooting as they have been 
influenced by the process for longer. Where possible, we have 
avoided short period systems as tidal forces could have altered 
the stellar structure and evolution. 


2 Select model sequences have been constructed using 1999 mesh 
points, and separately, twice the number of timesteps. The evolutionary 
tracks are indistinguishable from the ones presented. We are therefore 
satisfied that the computations are numerically converged. 
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System 

Period 

(d) 

Spectral 

type 

Mass 

(M 0 ) 

cr 

Radius 

(Rg) 

cr 

T e ff 

(K) 

cr 

logg 

cr 

log L/L 0 

cr 

V539 Ara 

3.17 

B3V 

6.240 

0.066 

4.516 

0.084 

18100 

500 

3.924 

0.016 

3.293 

0.051 



B4V 

5.314 

0.060 

3.428 

0.083 

17100 

500 

4.093 

0.021 

2.955 

0.055 

CVVel 

6.89 

B2.5V 

6.086 

0.044 

4.089 

0.036 

18100 

500 

3.999 

0.008 

3.207 

0.049 



B2.5V 

5.982 

0.035 

3.950 

0.036 

17900 

500 

4.022 

0.008 

3.158 

0.049 

V364 Lac 

7.35 

A4m 

2.333 

0.014 

3.309 

0.021 

8250 

150 

3.766 

0.005 

1.658 

0.032 



A3m 

2.295 

0.024 

2.986 

0.020 

8500 

150 

3.849 

0.006 

1.621 

0.031 

AIHya 

8.29 

F2m 

2.140 

0.038 

3.916 

0.031 

6700 

60 

3.583 

0.006 

1.443 

0.017 



F0V 

1.973 

0.036 

2.767 

0.019 

7100 

65 

3.849 

0.005 

1.242 

0.017 

El Cep 

8.44 

F3V 

1.772 

0.007 

2.897 

0.048 

6750 

100 

3.763 

0.014 

1.194 

0.030 



Flm 

1.680 

0.006 

2.330 

0.044 

6950 

100 

3.929 

0.016 

1.056 

0.030 

TZ For 

75.67 

G8III 

2.045 

0.055 

8.320 

0.120 

5000 

100 

2.908 

0.013 

1.589 

0.037 



F7IV 

1.945 

0.027 

3.965 

0.088 

6350 

100 

3.531 

0.018 

1.361 

0.033 

WX Cep 

3.38 

A5V 

2.533 

0.050 

3.996 

0.030 

8150 

250 

3.638 

0.005 

1.801 

0.054 



A2V 

2.324 

0.045 

2.712 

0.023 

8900 

250 

3.938 

0.006 

1.617 

0.050 

V1031 Ori 

3.41 

A6V 

2.468 

0.018 

4.323 

0.034 

7850 

500 

3.559 

0.007 

1.804 

0.112 



A3V 

2.281 

0.016 

2.978 

0.064 

8400 

500 

3.848 

0.019 

1.598 

0.105 

SZ Cen 

4.11 

A7V 

2.311 

0.026 

4.556 

0.032 

8100 

300 

3.485 

0.006 

1.904 

0.065 



A7V 

2.272 

0.021 

3.626 

0.026 

8380 

300 

3.676 

0.006 

1.765 

0.062 

AY Cam 

2.73 

A0V 

1.905 

0.040 

2.772 

0.020 

7250 

100 

3.832 

0.004 

1.280 

0.025 



F0V 

1.707 

0.036 

2.026 

0.017 

7395 

100 

4.058 

0.005 

1.042 

0.025 

AQSer 

1.69 

F5 

1.346 

0.024 

2.281 

0.014 

6430 

100 

3.850 

0.009 

0.901 

0.027 



F6 

1.417 

0.022 

2.451 

0.027 

6340 

100 

3.810 

0.012 

0.939 

0.042 


Table 1. Parameters of the eclipsing binary systems used in this study. All the data comes from the compilation of iTorres et al.l 

(f20Tol) . 


For each system, we run a grid of models across 4 metallic- 
ities (namely Z=0.01, 0.02, 0.03, 0.04) and with overshooting 
parameters <5 0V from 0 to 0.30 in steps of 0.03. The helium con¬ 
tent is presumed to vary as Y=0.25 + 15Z. To assess the quality 
of fit, for each timestep in the model sequence we calculate a 
goodness-of-fit via the formula: 



[Xj~ ft] 2 ' 

2of , 


( 2 ) 


where x,- is the relevant parameter from the stellar model (R, or 
T e fj) for each star in the system, and /./, and cr,- are the observed 
quantity and its error bar respectively. Note this means we are 
fitting both the primary and secondary simultaneously at each 
timestep. If our model fits perfectly we obtain P=l. We specif¬ 
ically use only R and T e g in this fit (and not surface gravity or 
luminosity) as these are the directly measured quantities. For the 
purposes of this study, we have not investigated the effects of 
uncertainties in the mass determinations. 


3.1. El Cep 

A best fit for El Cep is obtained with our Z= 0.02 model with 
d ov = 0.15 model, where we obtain P = 0.9066. In Fig. 2, we 
show the evolutionary tracks for both components of this system 
in the HR diagram, together with the observed characteristics of 
the system. Our best fit timestep is displayed by the crosses in 
the figure. 

For this system, we have also computed models for metallic- 
ities of Z=0.015 and Z=0.025 in order to get some idea of how 
sensitive d ov is to changes in metallicity. The best fit probabil¬ 
ity falls off rapidly as metallicity of the models is decreased. At 
Z=0.015, P = 3.292 x 10~ 2 , whereas the decline is much slower 
at higher metallicites, falling from P = 0.6880 at Z=0.025 to 
P = 0.3929 at Z=0.03. The reason for this is shown in Fig. 3. 
Overshooting has two effects: primarily it extends the length of 
the main sequence, but it also makes the star slightly more lu¬ 
minous at a given temperature. Lower metallicities favour lower 



Fig. 3. Evolutionary tracks in the radius-T e ff plane for the pri¬ 
mary star of El Cep computed with different metallicities: 
Z=0.015 (light grey), 0.02 (medium grey) and 0.025 (dark grey). 
Crosses represent the best fit models for each track and the er- 
rorbars denote the observed properties of the system. 


values of d ov because these models are already hotter and more 
luminous. The hook at the end of the main sequence occurs 
at higher temperatures at lower metallicity and for sufficiently 
small d ov , the observed temperatures cannot be reached while 
the star is on the main sequence. 

As an additional test for this system, we also compute mod¬ 
els where the initial helium content is varied by +0.05 (with a 
corresponding variation in the hydrogen abundance - i.e. the 
metal abundance, Z, is held constant). By reducing (increasing) 
the helium abundance, a model with the same Z and b ov be¬ 
comes larger (smaller) at a given temperature but the effect is 
quite small. This means that models with lower helium abun¬ 
dance tend to require smaller values of d ov . For Z=0.02, a model 


3 
























R. J. Stancliffe et al.: CUSP - calibrating convective overshooting 





Teff* 



Fig. 2. Left column: Evolutionary tracks in the radius-T e ff plane for the systems El Cep (top), AI Hya (middle) and WX Cep 
(bottom). The primary is denoted by the dark grey line and the secondary by the light grey line. Crosses represent the best fit models 
for each track and the errorbars denote the observed system. Right column: P-values for the systems as a function of age for the 
systems El Cep (top), AI Hya (middle) and WX Cep (bottom). The primary is denoted by the dark grey line and the secondary is 
denoted by the light grey line. The black line is the P-value for the system as a whole. 


with Y = 0.275 gives a best fit at 6 ov = 0.12 with P = 0.8595, 3.1.1. Comparison of evolution codes 

whereas one with Y = 0.285 gives <5 ov = 0.18 and P = 0.9092. 

Note that this is a marginally better fit than our standard case. 

In addition to our stars models, we have also computed a set 
of models using the stellar evolution code mesa j Paxton et all 
1201 lh . Following the prescription of lHerwiS (1200C ), convective 
overshooting in this code is implemented by means of a dif- 
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fusive exponential formalism whereby the extent of mixing is 
computed via the equation 


Dow = D 0 exp 



(3) 


where Do is the MLT diffusion coefficient inside the convec¬ 
tive region, z. is the distance from the convective boundary, H p 
is the pressure scale height at the convective boundary and / 
is a dimensionless free p arameter. Motivated by the work from 
IVandenBerg et al.l (2006), the overshooting parameter follows a 
ramp equation: 



1 - COS 71 


M, - M„ 


(4) 



Teff/K 


where /o is a constant, M* is the current stellar mass of the 
model, M m i„ is the stellar mass below which overshoot mix¬ 
ing does not occur, and M max th e stellar mass above which 
/ = fn. IVandenBerg et aD 1 2006 ) gives M m ; n = 1.1 M Q and 
M lnax = 1.8 M 0 . Previously, a parameter of/o = 0.014 was found 
to reproduce the width of the main sequence (for more details, 
see lHerwig||2000l) . 

The mesa model grid is run over the same metallicity range 
as the stars grid (namely Z=0.01, 0.02 and 0.03) and the free 
parameter fo is varied from 0 to 0.05 in steps of 0.005. 

The best fit model for El Cep as computed with mesa is for 
Z = 0.02 and fo = 0.04, with P = 0.399. As with the stars 
models, the quality of fit for the models falls off rapidly as the 
metallicity is changed such that no good fit is obtained at ei¬ 
ther Z = 0.01 or 0.03. In Fig. [4] we show the evolutionary 
tracks for the best fit models of both mesa and stars. The inter¬ 
nal structure of both stars for the best fitting models is displayed 
in Fig.0 The mesa model has a larger overshooting region than 
the stars model which is entirely consistent with the fact that 
the mesa evolutionary track has an end to the main sequence 
that is redward and brighter than the end of the main sequence 
in the stars track. If we express the extent of the overshooting 
region as a fraction of the pressure scale height at the convec¬ 
tive boundary, the primary and secondary for the stars model 
have fa = 0.251 and 0.244, respectively. For the mesa models, 
we obtain fa = 0.350 and 0.306 respectively. The stars hy¬ 
drogen profile shows much smoother features at the edge of the 
zone from which the fully mixed region has retreated. This is not 
due to any physical cause, but is a computational artefact. The 
stars code employs a non-Fagrangian mesh which is prone to 
numerical diffusion which acts to smooth out any composition 
discontinuities (for further details, see IStancliffejl2006l) . Above 
this region, the two models show very similar structure. 


Fig. 4. Evolutionary tracks in the radius-T e ff plane for El Cep. 
Solid lines denote models computed with stars, while dashed 
lines denote models computed with mesa. Dark grey denotes the 
primary and light grey the secondary. The best fit stars model 
is denoted with crosses, while the best fit mesa model is denoted 
with stars. The black errorbars denote the observed properties of 
the system. 


3.3. SZ Cen 

For SZ Cen, we obtain our best fit for Z=0.01 and <5 OV =0.12, with 
P= 0.7932. This places the primary beyond the end of the main 
sequence. iPolset alJ (11997) give a best fit metallicity between 
0.018 and 0.024, while IClared (120071) obtains a best fit at Z=0.02 
and with 0.1 pressure scale heights of overshooting. However, 
it should be noted that iTorres et al.l ( 2010 ) revised the param¬ 
eters for this system, compared to the values used in both of 
the above works. They adopted a temperature based on the beta 
index, rather than on an assumed spectral type coupled with a 
colour-T e ff calibration (G. Torres, private communication). The 
result of this is that the temperatures given bv lTorres et al.l d2010i) 
are somewhat hotter than those used in earlier works. If we adopt 
the parameters for SZ Cen used by Iciaretl ( 20071) . we obtain a 
best fit at Z=0.03 for 6 ow = 0.03 with P = 0.0121. This does not 
represent a good fit. On the other hand. [Ribas et al. ( 2000 ) obtain 
a best fit to this system with Z=0.007, Y=0.2 and 0.1-0.2 pres¬ 
sure scale heights of overshooting. Aside from the low helium 
value, our result is consistent with these parameters. 

3.4. AY Cam , CV Vel and V539 Ara 


3.2. VI031 Ori 

For V1031 Ori, our best fit system has a metallicity of Z=0.03. 
Both c) ov = 0.21 and 0.24 give the same P-value, namely 0.9297 
and it is p resumed that the true best fit lies between these two val¬ 
ues. |ciaret| (|2003) gives a best fit for Z=0.015 and 0.15 pressure 
scale heights of overshooting. iRibas et all (l2000h find similar 
parameters, with Z=0.016, Y=0.25 and over 0.2 pressure scale 
heights of overshooting. [Pols et all (119971) . on the other hand, 
prefer a model with overshooting and a metallicity of 0.023, 
though their model without overshooting and with a metallic¬ 
ity of Z=0.029 is also a good fit. It is difficult to reconcile these 
different models for the same system. 


These three systems all give excellent fits. For AY Cam, we 
obtain a best fit for Z=0.02 and <5 ov = 0.09, with a P-value of 
0.9267. We have been unable to find any other attempts to fit 
this system in the literature. CV Vel and V539 Ara are consid¬ 
erably more massive than the other binaries in our sample and 
both have primary masses around 6 M 0 . For both systems, we 
find a best fit at Z=0.02, though the S m values are different for 
both. For CV Vel, we prefer <5 nv = 0.09, while for V539 Ara d ov 
= 0.15 gives the best fit. [Pols et al.l (1997b also attempted to fit 
both these systems. For CV Vel, they give a slight preference 
to models without overshooting and suggest the best fit metal¬ 
licity is Z=0.016. For V539 Ara, their overshooting models are 
preferred and Z=0.016 is again the best fit metallicity. 
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Fig. 5. Hydrogen profiles as a function of mass for the best fit 
models for El Cep. The upper panel shows the primary star and 
the lower panel the secondary star. Dark grey lines represent the 
stars model, while the mesa model is in light grey. The star rep¬ 
resents the edge of the (Schwarzchild) convective core and the 
square the edge of the overshooting region. 



Teff« 


Fig. 6. Evolutionary tracks in the radius-T e jf plane for the com¬ 
ponents of SZ Cen. The primary is denoted by the dark grey line 
and the secondary by the light grey line. Crosses represent the 
best fit models for each track and the errorbars denote the ob¬ 
served system. 


Primary Secondary 

Mass (Mo) 3.86±0.15 2.95±0.09 

Radius (R 0 ) 32.9±1.5 3.0±0.2 

Te ff (K) 4 850±100 11000±500 


Table 2. Parameters of the system V2291 Oph. These data are 
taken from iSchroder et all (1 997 ). 


3.5. AI Hya, TZ For and V2291 Oph 

Though not included in the Torres et al. list, we have also at¬ 
tempted to model the system V2291 Oph. This C Auriga system 
was used by ISchroder et al. ( 1997 1 for their overshooting cal¬ 
ibration and it was a key system in their determination of the 
overshooting parameter. The parameters of this system are listed 
in Table [2] 

We are able to obtain reasonable, if not excellent, fits to the 
systems Al Hya, TZ For and V2291 Oph. These latter two sys¬ 
tems are significantly evolved. ISchroder et al] (119971) obtained 
6 m - 0 .12, for a model with X=0.70, Y=0.28 and Z=0.02. IClare1 
(2009^ chose a higher metallicity based on chemical analysis by 
iMarshalll d 19961) . preferring Z=0.03, and finds a best fit when 
0.2 pressure scale heights of overshooting are used. Our best fit 
is for Z = 0.03 and <5 0v = 0.24, with P = 0.4777. This places the 
primary on the first ascent of the giant branch. However, we also 
obtain a fit of P = 0.2240 for Z = 0.02 and 5 ov = 0.15 which 
would place the primary in the core helium burning phase, in 
agreement with the Pols et al. solution. 

TZ For is another system in which the primary is highly 
evolved. For this system, we obtain a best fit for Z = 0.03 and 
b ov =0.24, with P = 0.7036. This places the primary on the first 
ascent of the red giant branch. The P-values for neighbouring 
6 0v and Z-values in our grid are all less than 0.1, so this solu¬ 
tion is strongly favoured. Previous attempts to obtain a solution 
for this system have also had issues. IClareti d2007l) suggests us¬ 
ing 0.6 pressure scale heights of overshooting, but notes that this 
seems_to be too large for stars of masses comparable to TZ For. 
[Pols et all (119971) find slight evidence in favour of models with 
overshooting, but their solution is also less than ideal (see the 
right-hand panel of their figure 11). 

Unlike the previous two systems, Al Hya has both its com¬ 
ponents on the main sequence. For this system, we find a best fit 
for Z=0.04 and b m , =0.15. Fits from neighbouring metallicities 
are notably worse, as are neighbouring 6 0v values. Our fit is not 
perfect mostly because we do not find a good fit to the secondary, 
as shown in the middle-right panel of Figure [2] Our solution is 
in good agreement with iPols et all fl99 71). who found clear ev¬ 
idence for a superior solution when overshooting was included, 
but noted that the best fit metallicity lay outside the metallicity 
range (then' maximum Z was 0.033) of their isochrones. IClareti 
(12007! ) attempted to fit this system with models at Z=0.02 and 
was unable to obtain a satisfactory solution. 

3.6. V364 Lac, I NX Cep, AQ Ser 

For three of our systems, we are unable to obtain reasonable fits 
at all. For V364 Fac, WX Cep and AQ Ser, the best-fit P-value 
is below 0.1. 

Of the three systems, WX Cep is the least worst fitting sys¬ 
tem. We find a best fit value of P = 0.025 at Z=0.02 and b ov 
= 0.12. The main difficulty with this system is that the best 
fits to the components do not occur at the same age. The fit to 
the primary alone is at a maximum of about 0.75 at an age of 
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System 

Best fit Z 

Best fit (5 0v 

P 

V364 Lac 

0.02 

0.09 

6.914 x 10- 4 

AI Hya 

0.04 

0.15 

0.4922 

El Cep 

0.02 

0.15 

0.9066 

TZ For 

0.03 

0.24 

0.7036 

WX Cep 

0.02 

0.12 

2.549 x 10^ 2 

V1031 Ori 

0.03 

0.21(4) 

0.9297 

SZ Cen 

0.01 

0.12 

0.8767 

AY Cam 

0.02 

0.09 

0.9267 

AQ Ser 

0.02 

0.24 

8.501 x 10~ 55 

V2291 Oph 

0.03 

0.24 

0.4777 

CV Vel 

0.02 

0.09 

0.8632 

V539 Ara 

0.02 

0.15 

0.9267 



Fig. 7. Evolutionary tracks in the radius-T e ff plane for the com¬ 
ponents of AQ Ser, assuming a model with Z = 0.01 and <5 ov = 
0.18. The primary is denoted by the dark grey line and the sec¬ 
ondary by the light grey line. The errorbars denote the observed 
system. 

5.37 X 10 8 yr, but the fit to the secondary reaches a peak value 
of 0.92 at just 5.22 x 10 s yr. There is little overlap between the 
two fit distributions, resulting in a poor fit overall, as can be seen 
in Fig. [2] Our best-fit metallicity (and hence also the age of the 
system) agrees with that of Claret! (120071 ). who finds a best-fit so¬ 
lution for 0-0.2 pressure scale heights of overshooting. lPols et alJ 
( 1997 ) also favour a model with overshooting and a metallicity 
of Z=0.02. 

V364 Lac is a peculiar case. We are able to obtain almost 
perfect fits to the individual components of this system (with 
best-fit values of almost unity) at Z=0.02 and with d ov = 0.09. 
However, these solutions have very different ages. The primary 
is best fit at an age of 6.7 x 10 8 yr, while the secondary is best fit 
at an age of 5.9 x 10 8 yr. We are not aware of any other attempts 
to fit this system in the literature. 

AQ Ser has already been noted as a problematic system by 
iTorres et~aD (12014 ). These authors note that the more massive 
star appears systematically younger than its less massive coun¬ 
terpart. The P-value for our best fit to AQ Ser is so vanishingly 
small that it may well be considered zero. We can obtain a rea¬ 
sonable fit to the primary alone with Z=0.01 and h ov = 0.18 (see 
Fig. E), but it is clear from the evolutionary track that our sec¬ 
ondary is not massive enough. It is possible that a better fit to 
the system would be obtained if the masses of the two compo¬ 
nents were closer together, given how similar the radii and ef¬ 
fective temperatures of the two stars are. As a test of this, we 
ran a single calculation with Z=0.01 and c) ov = 0.18, in which 
we set the primary mass to be 1.395 M e and the secondary mass 
to be 1.370M o which are the minimum and maximum masses 
(respectively) allowed by the errors. In this case, we obtain a fit 
of P=0.101, which is a dramatic improvement over our original 
calculations. 

4. Discussion and conclusions 

We summarise the best fits to each of our systems in Table [2 
In Fig. [8] we have plotted the extent of overshooting as a func¬ 
tion of stellar mass, showing the overshooting both in terms of 
<5 ov and as a fraction of the pressure scale height. Note that we 
have not plotted TZ For and V2291 Oph on the pressure scale 


Table 3. Best fit properties for each of our systems using only T 
and R. 


height plot, as these systems are evolved beyond the main se¬ 
quence and we cannot define this quantity if there is no longer 
a convective core. For all of our systems, we require some de¬ 
gree of overshooting and in agreement with previous works the 
amount required is only moderate. The extent of overshooting 
ranges from 6 m = 0.06 to 0.24, though the highest values are 
found only in our most evolved and most problematic systems 
(AQ Ser, TZ For and V2291 Oph). Most systems have b ov be¬ 
tween 0.09 and 0.15. This is consistent with the 0.1-0.3 pres¬ 
sure scale heights of overshooting found by other authors. We 
see no reason to suggest that 6 ov is a function of mass (unlike 
iRibas et aO2000 ). though we note that the majority of our sys¬ 
tems are clustered around 2 M 0 , with just two systems at around 
6 M 0 . There is also no evidence for a metallicity dependence to 
the extent of overshooting, though this may simply be because 
we lack enough systems across a range of metallicities. While 
our three Z=0.03 systems do have higher 6 0v values, we stress 
that two of these are very poor fits hence we cannot provide a 
definite conclusion regarding a possible metallicity dependence. 

The extent of overshooting we find is consistent with values 
obtained from asteroseismic determinations. IXer tsj d2015l) finds 
values between 0-0.5 pressure scale heights, also with consider¬ 
able star-to-star variations. She also reports that there is no ob¬ 
vious relation between extent of overshooting and stellar mass. 
iNeiner et al.l (2012 ) find 0.3-0.35 pressure scale heights of over¬ 
shooting in their two late Be stars, HD 181231 and HD 175869. 
Again our results are consistent with these values. 

To complement our stars models, we have also calculated 
best fit models for each of our systems using mesa. The details 
of these fits are given in the appendix. We find general agreement 
between the two codes. Both predict the same metallicity for all 
but one of the systems (the exception is AI Hya). There is also 
agreement in the trend of the overshooting parameter required: 
for systems where stars requires a high 6 m , we also find that 
mesa requires a high fo. The extent of overshooting in the mesa 
models is about 0.15-0.4 pressure scale heights. Computing a P- 
weighted average of the overshooting parameter /o results in a 
value of 0.031. 

Ideally, one would like to be able to select a single value for 
the overshooting parameter that would work reasonably well for 
all systems. This is particularly relevant for future work in our 
project, where we aim to create grids of low-mass stellar mod¬ 
els for use in bonnsai. If we neglect the systems with P-values 
less than 0.1 and take an average of the overshooting param¬ 
eter for the remaining systems, in which we weight each sys¬ 
tem by its relative P-value, we obtain 6 ov = 0.156. Adopting this 
value, we would still obtain good fits for the systems CV Vel and 
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Fig. 8. Extent of overshooting as a function of mass. The upper 
panel displays the overshooting in terms of c)„ v , while the lower 
panel shows the same quantity in terms of a fraction of the pres¬ 
sure scale height. Symbols denote the best fit metallicity for the 
systems: Z=0.01 (triangles), Z=0.02 (squares), Z=0.03 (circles) 
and Z=0.04 (pentagons). 

V1031 Ori (which give P=0.7551 and 0.7243 respectively for<5 ov 
= 0.15). AY Cam and V2291 Oph would be moderate fits (P = 
0.2733 and 0.2240 respectively), though the best fit metallicity 
for the latter system would drop to Z = 0.02. For the systems SZ 
Cen and TZ For we are unable to obtain any fit at 6 ov =0.15. Of 
the nine systems we are able to fit when our choice of Z and 6 ov 
is free, we would only fit seven of these systems if we adopted 
our average value of d' ov = 0.156. 

It would be useful to be able to extend this work to a wider 
sample of binaries, particularly as many of our systems cluster 
around 2 M 0 . Additional data for stars of between 3 and 5 M 0 
(and above 6M Q ) would be particularly desirable. However, as 
can be seen from Fig. [7] few systems in the lTorres et ali d20ld) 
data set lie in this mass range even before we take into account 
our requirements that the systems should be fairly wide and con¬ 
tain components that are evolved away from the ZAMS. Nature, 
it seems, has not been kind to us in this regard. One can hope 
that future survey missions, such as Gaifl may help to fill this 
deficiency. 


System 

Best fit Z 

Best fit fo 

P 

V364 Lac 

0.02 

0.030 

2.244 x 10~ 3 

AIHya 

0.03 

0.030 

1.604 x 10- 3 

El Cep 

0.02 

0.040 

0.3987 

TZ For 

0.03 

0.050 

0.01072 

WX Cep 

0.02 

0.030 

0.05469 

V1031 Ori 

0.03 

0.045 

0.9474 

SZ Cen 

0.01 

0.025 

0.01779 

AY Cam 

0.02 

0.020 

0.5296 

AQ Ser 

- 

- 


V2291 Oph 

0.03 

0.045 

0.09719 

CV Vel 

0.02 

0.030 

0.9427 

V539 Ara 

0.02 

0.035 

0.9808 


Table A.l. Best fit properties for each of our systems using only 
T and R as determined using the mesa code. No reasonable fit 
could be found to AQ Ser (P < 10~ 96 ). 
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Appendix A: mesa model fits 

In addition to El Cep, we have also computed models for all 
our systems using mesa. The best fit parameters are given in 
Table lATI Fig. IA. II shows this data, together with the extent of 
the overshooting region as a function of the pressure scale height 
at the top of the convective region. Calculating an average /o 
weighted by the P-values, we obtain fo = 0.031. 
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Fig.A.l. Extent of overshooting as a function of mass for the 
mesa models. The upper panel displays the overshooting in terms 
of 6 m , while the lower panel shows the same quantity in terms 
of a fraction of the pressure scale height. Symbols denote the 
best fit metallicity for the systems: Z=0.01 (triangles), Z=0.02 
(squares), Z=0.03 (circles). 
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